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Abstract 



Using numerical ray tracing, the paper studies how the average distance modulus in 
an inhomogeneous universe differs from its homogeneous counterpart. The averaging is 
over all directions from a fixed observer not over all possible observers (cosmic), thus it 
is more directly applicable to our observations. Unlike previous studies, the averaging is 
exact, non-perturbative, and includes all possible non-linear effects. The inhomogeneous 
universes are represented by Sweese-cheese models containing random and simple cubic 
lattices of mass-compensated voids. The Earth observer is in the homogeneous cheese which 
has an Einstein - de Sitter metric. For the first time, the averaging is widened to include 
the supernovas inside the voids by assuming the probability for supernova emission from 
any comoving volume is proportional to the rest mass in it. For voids aligned in a certain 
direction, there is a cumulative gravitational lensing correction to the distance modulus that 
increases with redshift. That correction is present even for small voids and depends on the 
density contrast of the voids, not on their radius. Averaging over all directions destroys the 
cumulative correction even in a non-randomized simple cubic lattice of voids. Despite the 
well known argument for photon flux conservation, the average distance modulus correction 
at low redshifts is not zero due to the peculiar velocities. A formula for the maximum 
possible average correction as a function of redshift is derived and shown to be in excellent 
agreement with the numerical results. The formula applies to voids of any size that: (1) 
have approximately constant densities in their interior and walls, (2) are not in a deep 
nonlinear regime. The actual average correction calculated in random and simple cubic void 
lattices is severely damped below the predicted maximum. That is traced to cancelations 
between the corrections coming from the fronts and backs of different voids at the same 

viii 



redshift from the observer. The calculated correction at low redshifts allows one to readily 
predict the redshift at which the averaged fluctuation in the Hubble diagram is below a 
required precision and suggests a method to extract the background Hubble constant from 
low redshift data without the need to correct for peculiar velocities. 
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1 Introduction 



It is fair to say that the standard cosmological LCDM model is facing a phenomeno- 
logical crisis. The dark matter carrier has been evading direct detection for decades and 
the origin of dark energy remains a theoretical puzzle. The most natural candidate is the 
vacuum state energy, but the flat-space Quantum Field Theory is incapable of calculat- 
ing it, producing an estimate that is 10 120 times as large as the value suggested by the 
supernova observations. Presumably, that would be rectified in a fully quantized theory 
of gravition which unfortunately does not yet exist. Remaining within the standard Gen- 
eral Relativity, there are attempts to explain dark energy as an apparent quantity arising 
from the averaging procedure that maps the real lumpy universe to the idealized homo- 
geneous Friedman-Lemaitre-Robertson- Walker (FLRW) metric. Such approaches separate 
into three distinct classes. First, one could average the Einstein equations over spatial 
hypersurfaces and hope the analogue of the Friedman equation contains a significant effec- 
tive lambda term. Such an idea is conceptually problematic since inhomogeneous universes 
do not naturally pick up preffered spatial slices to average over, unlike the homogeneous 
FLRW models. It is shown in [l] that the result of such averaging depends on the arbitrary 
choice of slicing and have no coordinate independent meaning - even Minkowsky spacetime 
could produce an apparent acceleration. A cousin to the first_ap_proach is averaging over 
the null hypersurface of the past light cone of the observer, 
such averages have a physical meaning and are not simply coordinate quantities. In the 
second approach, the second order perturbations to the Einstein equations are calculated 
and treated as an effective energy momentum tensor term. It is pointed out in jl| that the 
effective tensor arising from perturbation theory is not gauge independent and therefore 
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3). It is not proven that 



cannot be used as an effective energy momentum tensor in the Einstein equation. 

The present paper belongs to the third approach which is completely coordinate inde- 
pendent since it calculates relations only between physically observable quantities such as 
redshift and luminosity distance (or its logarithmic measure - the distance modulus). Pa- 
pers of this type, for example j4j, demonstrate that the distance modulus - redshift Hubble 
diagram of type la supernovas can be reproduced without any dark energy for an observer 
occupying the center of a giant (~ Gpc) under-dense spherical bubble. That is possible 
since the central observer measures a local Hubble parameter that is larger than the global 
one. The corresponding shift on the Hubble diagram between the true background and the 
one assumed by the observer allows for fitting the supernova data. Bubbles of such a large 
scale have significant peculiar velocities away from their center. That creates a fine tunning 
problem - the observer has to be improbably close to the bubble center 5j to observe a 
CMB dipole in agreement with the measured v/c ~ 0.002 (Local Group velocity ~ 620 
km/s). This violates the Copernicus principle that we do not occupy a preffered position 
in the universe. 

Other papers of the same class studied the cumulative effect of a configuration of 
many bubbles/voids of smaller size, thus avoiding big peculiar velocities and significant 
anisotropics in CMB. That scenario can be conveniently simulated in the "Swiss-cheese" 
toy model jf], Q, 0] which has the virtue of being analytically solvable. It is constructed 
by removing spherical regions from a homogeneous FLRW background (the "cheese") and 
replacing them with inhomogenous density distributions with the same gravitating mass 



(mass-compensated voids) . It was used in 



9|, yj for an observer looking along the diameters 



of a string of aligned voids of radius 350 Mpc. A positive cumulative change in the dis- 
tance modulus was found which increased with redshift with respect to the homogeneous 
background. Although the size of the correction was significant at high redshift, it was 
not large enough to fit the supernova Hubble diagram and it substituted only partially for 



dark energy. A smaller correction due to aligned voids with a more mundane radius of 



23/i" 1 Mpc was calculated in 



10]. The result obtained in sj] was criticized in 
was demonstrated to stem from the cumulative weak lensing defocusing of the rays passing 
diametrically through the aligned voids. The same paper showed that the correction to the 
distance modulus goes to zero when averaged over randomized impact parameters of the 
rays entering the voids. This conclusion is in accord with an old argument [12j based on 



11] where it 



the gravitational lensing conservation of the total photon flux. The implicit geometrical 
assumptions in [ljj] have been challenged in more recent papers Q] (and references 
therein) and the present work will demonstrate they are violated close to the observer due 
to peculiar velocities modifications of the redshift surfaces. 

The studies mentioned above are defficient in several areas. Averaging over directions 



was not performed in 



llj were done with weak lensing theory 



9|. The calculations in 
neglecting the time dependence of the void density and possible nonlinear effects which 
become significant at low redshifts. The averaging assumed that the ray impact parameters 
had a uniform probability distribution over the cross-section of a void. That is increasingly 
true at high redshifts but doesn not hold in the local neighbourhood where the rays have 



11| averaged only over supernovas residing 



to converge on the observer. More importantly, 
in the homogeneous cheese. The same was done in a study that obtained an analytical 
estimate of the corrections to the luminosity distance, [15]. Such a bias is unjustified, 
first because supernovas in the voids produce much bigger corrections than the ones in the 
cheese 9] , and second the observed supernovas occur more frequently in denser regions like 
the void walls than in the cheese. 

So far, the question whether the direction-averaged distance modulus correction from 
a configuration of voids really vanishes has not been answered in the context of a calcu- 
lation that: (1) is exact, non-perturbative, an includes all possible non-linear effects; (2) 
is armed with a physically sensible averaging procedure free of ad hoc assumptions about 



121 ] ; (3) does not neglect the 



impact parameters [11] or random cancelations of corrections 
supernovas inside voids. 

The goal of the present paper is to rectify that situation. The final result of the cal- 
culations can not be guessed on the grounds of the previous studies since it may depend 
on non-linear effects and the choice of averaging procedure and supernova selection none 
of which was taken into account before. The calculations are performed in Swiss-cheese 
models with mass-compensated voids of two radii having observational support: 30 and 300 
Mpc. The cheese is spatially-flat matter-only Einstein - de Sitter (EdS). The observer is 
placed in the cheese since at present there is no indication that we occupy a void. The first 
Swiss-cheese model is set up in section 2. The observer shoots past-directed light rays in all 
directions. The light propagation geodesic equations and the luminosity distance tracing 
along each ray are discussed in section 3. The physically sensible averaging procedure in 
this paper assumes that the probability for supernova emission from a comoving volume 
is proportional to the rest mass in it. Averaging the distance modulus for a single void is 
discussed in section 4 and it carries the essential characteristics of the procedure for many 
voids. The main outcome of that section is a simple formula to estimate the maximal aver- 
age corrections to the distance modulus. Section 5 studies the cumulative correction due to 
gravitational lensing along a string of aligned voids. Numerical results from averaging over 
random and simple cubic void lattices are presented in section 6. The effect on the distance 
modulus produced by large voids of radius 300 Mpc, that leave a measurable imprint on 
CMB, is evaluated in Section 7. The summary and conclusions section discusses the answer 
to the main question addressed by this paper and the practical importance of the obtained 
low-redshift results for future surveys. 

The speed of light in this paper is c = 1, and times and distances are measured in 
mega parsecs (Mpc), occasionally giving time in mega years (Myr) for convenience. The 
gravitational constant G is kept explicit in all equations so that the reader can easily 



substitute with a favorite value. The usual geometrized units, G — 1, were used in the 
numerical calculations. 
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2 Swiss-cheese model with voids of radius r v = 30 Mpc 



2.1 Einstein-de Sitter as homogeneous cheese 

past while the structure 



161 ] on the development 



Every reasonable model of the universe is matter dominated in the 
is being formed. Adding a cosmological constant has little impact p 
up until recent times when it starts dominating the matter. That motivates choosing the 
spatially-flat matter-only Einstein-de Sitter (EdS) model for the homogeneous regions of 
the Swiss-cheese. The matter includes both visible and dark varieties. This is a convenient 
playground to explore the question whether voids in the matter distribution can mimic 
the effect of dark energy. The Hubble constant of the model is Hq = 70 km/s/Mpc. 
The big bang time is conventionally chosen as = 0, the current age of the model is 
to = 2/(3Hq) = 2857 Mpc (9312 Myr) and the scale factor as a function of cosmic time 
is a(t) = (t/to) 2 ^- The matter density today is p(io) = 3Hq/(8itG) and as a function of 
time is: 

m = E H = ^- (2.i) 

2.2 Lemaitre model of a spherical void 

The voids in the Swiss-cheese are modeled with the Lemaitre metric 6_|, |7J also known as 
" Lemaitre- Tolman-Bondi solution". It describes a spherically symmetric spacetime filled 
with an irrotational pressureless ideal fluid (matter, dust). The matter particles are in a free 
fall under their own gravity tracing geodesies. The zero rotation of the geodesic congruence 
means the geodesies are hypersurface-orthogonal and the corresponding family of spatial 
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hypersurfaces define a convenient foliation and coordinate system on the spacetime. The 
resulting coordinates = (t, r, 8, 0) are matter-comoving and synchronous. The matter 
energy-momentum tensor in these coordinates is diagonal: T^ v = diag(p(r,t), 0, 0, 0). The 
metric is given by: 

ds 2 = -dt 2 + — dr 2 + R(r, tf(d9 2 + sin 2 9 d<j) 2 ) , (2.2) 

1 + 2E[r) 

where prime denotes a derivative with respect to the radial coordinate r. The arbitrary 
integration function E(r) results from integrating the Gf r = Einstein equation. The areal 
radius R(r,t) determines the area of a sphere of radius r and E(r) determines the local 
3-curvature of the spatial slices. The time coordinate t measures the proper time of the 
comoving matter. Integrating once the G rr = Einstein equation leads to the following 
evolution equation: 

where dot denotes a derivative with respect to t and M(r) is another arbitrary integration 
function. The above equation alludes to conservation of kinetic plus gravitational energy 
in Newtonian mechanics where the "energy function" E(r) plays the role of total energy 
per unit mass. The rest of the Einstein equations connect the "mass function" M(r) to the 
comoving matter density p: 

M'(r) F 
£ ) = a n2 ( \L ( a > M ( r ) = / 47rR (^ VPir, t)df . (2.4) 

The integral can be evaluated for any t and will produce the same M(r). Note, M(r) differs 
from the comoving rest mass which has an extra factor of (1 + E{r))~^-I 2 in the integral. 



This is an example of a "relativistic mass defect' 



3, [8] . The parametric solution of (12. 3ft 
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for the E > case of interest is: 



R(r,t) 



GM(r) 
2E(r) 



(cosh rj — 1) 



t-t B (r) 



(2E(r)) 3 / 2 



GM(r) 



(sinh rj — rj) 



(2.5) 



where the "bang time" tg(r) is another arbitrary function and rj > is a parameter. It 
is natural for inhomogeneous models to have different places with big bang happening at 
different times. 

The metric ( 12.2ft contains as a subclass all FLRW models which are obtained by setting 
R = a(t) r and E = —kr 2 , k = const. In particular, the cheese regions of the present model 
are spatially-flat and have E = 0. The voids have to be matched to the homogeneous EdS 
without tearing of the metric. The appropriate junction conditions are: 



where r v is the radius at which the void merges with a homogeneous region. The first two 
conditions express an obvious continuation to the EdS values. The third condition means 
the void is " mass-compensated" : the gravitating mass inside it matches the mass in EdS 
within the same radius r v . The conditions above guarantee that the metric outside radius 
r v remains EdS exactly - it does not " feel" the presence of the void and the metrics inside 
different voids do not interfere with each other. Of course, voids in the real world are 
not restricted within some radius, they interfere and even merge, their walls accrete mass 
from outside and regions where the metric stays exactly homogeneous do not exist. The 
inconvenience is that a real world simulation requires computing of the global metric (or 
its Newtonian approximation) encompassing the whole spacetime. The mass-compensated 
voids are a toy-model which is easier to compute because the metrics inside different voids 
are exact copies of each other. 

The Lemaitre metric is specified by choosing three functions, say R(r,ti), E(r), and 




E(r v ) = 0, 



t B (r v ) = 0, 



M(r v ) = M EdS (r v ) , 



(2.6) 
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M(r), where t{ is some initial time and the last two functions determine R{r,t) via (12. 3p . 
The radial coordinate r is just a label for the spherical matter shells and can be chosen 
at will (gauge freedom). That means the system is specified by only two functions plus a 
gauge choice for r, effected by choosing R(r, tj). Various combinations of two functions are 



discussed in 



17| 



The r gauge in the present paper is fixed by choosing 

R(r,t ) = r. (2.7) 



It was demonstrated in 18[ that the Lemaitre model can be written from synchronous to 



perturbed FLRW coordinates (known as "Newtonian gauge"): 

ds 2 = -(1 + 2ip)dt 2 + a(t) 2 (l - 2tp)(dx 2 + dy 2 + dz 2 ) , (2.8) 

as long as the peculiar velocities remain small (which is demonstrated later). The void 
model needs initial conditions in the synchronous coordinates (I2.2p but the astronomical 
data uses the Newtonian coordinates (12.81) . This dichotomy is resolved by the choice in 



18| 



(12. 7p which makes the two coordinate systems coincide, approximately, at time to 
That allows for using Newtonian notions in synchronous coordinates around time to and 
feeding astronomical data of distances, densities and velocities directly to model. Other 
models in the literature fix the synchronous r gauge at some initial time much earlier 
than tQ. By the time to, their synchronous coordinates evolved significantly away from the 
Newtonian ones which prevents direct Newtonian interpretations and may lead to illusory 
coordinate effects. For example, light rays that look straight in Newtonian gauge often 
appear "curved" in such synchronous coordinates. 

The present paper uses p(r, tg) and tg(r) as the two functions specifying the void metric. 
There is not a standardized definition of what constitutes a void in astronomy. The so called 
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supervoids are defined as regions free of rich clusters and in the Milky Way proximity have 



an average radius of about r v ~ 46 h 1 Mpc ~ 65 Mpc (for Hq = 70 km/ s/Mpc) 



19|]. 



Computer void-finding algorithms usually define a void nucleus to be completely free of 



20] 



galaxies leading to breaking down the supervoids into smaller ones of r v ~ 20 Mpc 
Different algorithms often find voids of different sizes, shapes and numbers for the same 
region of space 21J. A recent paper examined the SDSS 5th release data 22J using an 



algorithm that finds voids closest to visual inspection and found an average r v ~ 36 Mpc 



22] . It is obvious there is not an universal agreement upon the average void size and the 
current paper adopts an intermediate value of r v = 30 Mpc. The conclusions reached at 
the end scale with r v as long as the voids are not significantly nonlinear. The current 
galaxy density inside voids is very low, below 10% of the average. Density reconstruction 
of our neighbourhood using the peculiar velocity field shows that the total density (dark 
and visible matter) inside voids dips to below 0.40 p(io) [23|- The bias between visible and 
dark matter is still an open question but the present paper will assume that light is a good 
tracer for matter and set the total matter density inside the void to Pi n (to) = 0.10 p(to). 
Observationally, most of the visible matter is gathered in the void walls and it will be 
assumed the same applies to the dark matter distribution as well. The wall thickness is set 
as ~ 5 Mpc. The so designed matter density of the void model is 



Ax + A 2 tanh[2(r - 25)] - A3 tanh[3(r - 29)] , r < 30 
p(r,t ) = p(t ) { (2.9) 

, r > 30, 



and is shown on Fig. 12.11 The values of the constants to six significant figures are 
(At, A 2 , A3) = (0.548006, 1.25446, 0.806455). They were determined by requiring den- 
sity continuity and satisfaction of the compensated mass condition in (12.61) . The mass 
function is obtained by evaluating the integral in (12 .4p at time tQ. 
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P (r, t ) I p (t ) 
2.5 : 
2 . 
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1 . 
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10 15 20 25 30 35 



Figure 2.1: Current void density 
p(r, to) divided by the homogeneous 
r ' Mpc density p(to). 



The Lemaitre metric is fixed comp 
be the peculiar velocity at time 



etely by specifying one more function. This could 
171 ]. but it is not clear what a typical velocity profile 
looks like. The bang time t B {r) turns out to be a better candidate. After specifying it, 
one needs to solve for E{r) the two equations in (12.51) evaluated at time to- Eliminating r] 
gives 

((1 + A x f - l) 1 / 2 - arccosh(l + A x) = x 3/2 (t - t B (r)) , (2.10) 

where A (r) = R(r,t )/(GM) 1 / 3 and x(r) = 2i?/(GM) 2 / 3 . The condition for equation 
(I2.10p to have a solution x(r) at a given r is derived in [171 ] : 



t 



-s( r ) > tcritir) = 



. V2A (rf/ 2 
= t § • 



(2.11^ 



Equation fl 2 . 1 [) can be solved numerically for x(r) which in turn gives E(r). All the models 
with the given density (12.91) are enumerated by the possible functions ig(r) > t cr ^{r). The 
present model sets 

t B (r) = 0, (2.12) 
which satisfies the junction condition (I2.6P and the solvability condition (12. lip . 



Choosing dt B /dr = excites only the growing density perturbation mode 



24 



25] 



which is favored by linear perturbation theory in cosmology and increases like 6 = Sp/p 
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oc a(t) in EdS [26;]. The decaying mode 5 oc H(t) 



a/ a [27| decreases like 



oc t" 1 in EdS. It is neglected in cosmology assuming it had enough time to decay to 
insignificant levels. Choosing dtj^ jdr 7^ excites a mixture of growing and decaying modes. 
A common example are Lemaitre models that set the peculiar velocity at initial time 
to zero which in effect creates a sum of growing and decaying modes with amplitudes 



Ugrow decay (U) = 3/2 

The model selected by (12.121) is a pure growing mode. Knowing E(r), the metric 
function R(r, t) can be obtained by either using the analytic solution (12.51) or by numerically 
integrating (I2.3p . The second method is advantageous in Mathematica since it produces 
an interpolation function which is fast to call and whose numeric derivatives are easy to 
calculate. 

2.3 Void model properties: density, shell crossing, last scattering, 
peculiar velocity 

The comoving matter density of the void is shown on Fig. 12.21 and behaves as expected. 
The ratio of the void wall density to the homogeneous density is increasing at all times. 
The wall eventually collapses, reaching infinite density at shell crossing which happens at 
t « 4500 Mpc (14680 Myr). 




p(r,t) / p(t) 




r, Mpc 



Figure 2.2: Density ratio p(r, t) / pit) at 
times (500, 1000, 2000, t , 3500, 3900) 
Mpc (from bottom to top). 
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Figure 2.3: Density ratio p/p at time 
r, Mpc Q £ j ag ^ sca ^ erm g 



The time of last scattering in this model, defined as the time at which the EdS scale 
factor is a = 1/1100, is t^g = 0.0783 Mpc (0.255 Myr). The density ratio at last scat- 
tering is shown on Fig. 12.31 The density contrast of the perturbation is of the expected 



order: 5 = p/p — 1 



10 



-3 



This 5 applies to the dominating dark matter; the baryonic 



matter perturbations, suggested by the CMB data, are still 5i ar ~ 10 5 at the time of last 
scattering since the baryons had just decoupled from the photons. 
v(r, t) 




Figure 2.4: Peculiar velocity profiles at 
r M c times t=(t LS , 10, 100, 1000, to, 4500) 
35 ' Mpc (from bottom to top). 



The peculiar velocity is defined as the matter velocity with respect to the Newtonian 
coordinates ( 12. 81) . This concept does not apply to synchronous coordinates at all but can 



be calculated from quantities in them, 



18|: 



v(r,t) = R(r,t) - R(r,t) 



a(t) 



v < 1 



(2.13) 
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28|, 



Fig. 12.41 shows a few velocity profiles at different times. The peculiar velocity is everywhere 
positive since in Newtonian coordinates (12. 8p the matter underdensity inside the void acts 
like a repelling gravitational source which accelerates the matter at bigger radius outwards. 
Even if the initial velocity was negative (towards the center), it will be reversed at later 
times outwards. For an underdense void interior, a negative initial velocity there is possible 
only if the model contains a decaying perturbation mode. The current model has a pure 
growing mode and the peculiar velocity is always positive. 

The velocity peak in the profiles can be estimated with standard linear perturbation 
theory. The linearized mass continuity equation in Newtonian coordinates gives 
V -v — —ad5/dt = — a 5 dln5/dlna. If the void interior is homogeneous, the velocity 
there is linear (see Fig. I2.4D : v(r,t) = ft(t)r, where (r, t) are Newtonian coordinates (12. 8ft 
comoving with the averaged matter, not the synchronous ones. For the growing mode, 
5 oc a in the EdS. Plugging back in the continuity equation obtains 3k = —a 5 which leads 
to 

v = --ar5 = --Har6. (2.14) 
3 3 K J 

Fig. 12.41 shows that the velocity peaks at synchronous r ~ 25 Mpc. Taking that value 
for the Newtonian r and noting that today 5(to) = —0.9 in the void interior one gets 
v max{to) ~ 1.75 x 10 ^ when the actual value from Fig. 12.41 is v m ax 

(t ) ~ 2.26 x 10 _ <\ 

The agreement is not bad if one notes that the perturbation is already in nonlinear regime 
at time to- 

The shape of the velocity profiles is explained by the velocity asymptotic behavior close 
to the big bang when t — tg(r) — > and 77 — >■ 0. One can expand the second equation in 
( 12.51) and then invert the series to obtain 



V 



61/3b -^ 3+ ^ b5+0(b7) - bm » - ^ ( w) V3 ■ (215) 
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Substituting that in the first equation gives 



„, , GM/3 2 /3 3 4 / 3 
2E V 2 1 / 3 2 2 / 3 10 



5 4 - 



27 
1400 



E 6 + 0{B h 



(2.16) 



Inserting the above equation into the peculiar velocity definition ( 12. 131) and using a/ a 
2/(3*) for EdS obtains 



v(r,t) = 2E 



t-t B y/ 3 ({A/^/Hb , (3/4)V3(i + tfl / t ) 9(2 + * s /*) 



GM 



5 2 * 



700 



5 2 + 



(2.17) 

l lH B /t~t-^ 



26| . As mentioned 



The first term in the above expansion has a time dependence (* — *g)" 
for small *# <C *. One can identify it with the decaying density mode, 
before, that mode vanishes for models with t B = const (=0 by (I2.6P ). In this case, for 
times not too far in the future: * <C GM / {2E)^I^ (ensuring B < 1) , the second term 
in the above expansion dominates and the velocity profile evolves in time as ex 

*V3. This 



corresponds to the velocity of the growing density mode, [26J|. The constant shape of the 
velocity profile is set by E(r)/(GM(r))^l^ which is typically triangular since that quantity 
is zero at r = and r = r v and positive in between for E > models. 
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3 Ray tracing 



3.1 Redshift 

Light rays are null geodesies. The 4-momentum of a photon along such a geodesic is 
K.V = ko (dt/dX, dr/dX, d8/dX, d<p/dX) were kg is a constant depending on the choice of 
the affine parameter A. The redshift z of a photon emitted by a supernova and observed 
on Earth satisfies 

1 , „ Kbs 

= {U ■ K) em = {dt/dX) em 

^em ^obs 

(U ■ K) ohs (dt/dX) 

obs 



where is the 4- velocity of emitter and observer in the comoving coordinates (I2.2p and dot 



denotes scalar product. Both the supernova and the Earth observer are assumed comoving 
with the matter therefore = (1, 0, 0, 0) which leads to the last equality. The affine 
parameter is uniquely defined up to scaling and shifting (affine transformations). One can 
use those to chose a convenient parameter that satisfies 

X obs = 0, (dt/dX) obs = 1 , (3.2) 

where obs means evaluated at the Earth observer. Substituting that in ( 13. ip gives the 
connection between time and redshift along a light ray: 

^ = »M+1. (3-3) 
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3.2 Local and global spatial coordinates for ray tracing 

The affine parameter A is chosen to increase with the time t along the ray. The ray's 
geodesic equations are best integrated by following the ray back in time (decreasing A) 
from the Earth observer (A = 0) to a supernova (A < 0). Integration forward in time would 
require knowing the precise ray direction at the supernova required to hit the Earth and 
that is hard and impractical to calculate. Note that tracing a ray by decreasing A does not 
flip the sign of any derivative with respect to A! 

The fastest method of ray integration inside a void is to use local polar coordinates 
(r, 9) in the two-dimensional plane containing the ray and the void center, Fig. 13.11 Due 
to the symmetric metric (12. 2p . a light ray always remains in that plane which can be 
taken to be = const by properly orienting the spatial coordinate axes. This reduces the 
number of independent variables by one and simplifies all equations. From void to void, 
the rays are traced in global spatial coordinates comoving with the homogeneous matter. 
In the homogeneous cheese, these coordinates are the usual FLRW Cartesian 3-vectors 
x = (X, Y, Z). Inside a void (for r < r v ), the global coordinates are Cartesian-like and are 
connected to the local (r, 9) by the familiar relation: 

x = x c + r cos 9 u + r sin 9 v , (3.4) 

where x c contains the global spatial coordinates of the void center C and the unit 3-vectors 
{u, v} form a basis in the two-dimensional plane of the ray, Fig. 13.11 The above equation 
does not carry the usual physical sense of conversion from polar to Cartesian coordinates - 
the spatial metric inside a void is not really Cartesian if written in the global coordinates 
(13 .4p . Despite that, the coordinate transformation (13 .4p is mathematically permissible and 
one is free to use it. 
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Figure 3.1: Local coordinates in the 
ray plane: C - void center, F - point 
of ray exit, {u,v} - basis unit vectors. 



The u vector is chosen to point from the void center to the point F at which the ray 
leaves the void (enters the void when propagating back in time). The angle 9 is measured 
from the u-axis. The v vector is chosen orthogonal to u pointing in the half-plane occupied 
by the ray. The global components of {u, v}, needed to calculate the components of x, are 
determined geometrically from the global coordinates of the void center, x c , the point F, 
and the ray direction at F. 

A ray trajectory consists of segments inside the homogeneous cheese and segments 
inside voids which are glued together in the global coordinates ( 13. 4p . The final conditions 
of a cheese- segment determine the initial conditions for integration inside the following 
void-segment. 
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3.3 Void-segments 



The geodesic equations in the plane of the ray (0 = const) are 



d 2 t R'R' fdr\ 2 ■ / d6^ 2 



+ rr^br + RR l^) = ( 3 - 5 ) 



dA 2 1 + 2E V^A/ V rfA 
d 2 r R'dtdr fR" E' \ ( dr\ 2 . R f d6\ 2 

^ + 2-—— + 2——— = — ( R 2 —) = 0, (3.7) 



d\ 2 RdXdX RdXdX dX \ dX y 
where A is the affine parameter along the geodesic. The null condition is the first integral 

«' 2 '*)V*v wa\ (3 .8) 



1 + 2E \dXJ \dXJ \dX 

To decrease the computational time, the total differential order of the system can be reduced 
by using first integrals: 

TX = + -(*/*>') -j^ 0.9) 



2+1 



dt 
dX 

dr_ V1 + 2E 

dX ~ W 
d9 L 

dX ~ 'R 2 ' 



z + l) 2 - (L/R) 2 



where the last equation is an integration of (13.71) and L is the integration constant analogous 
to angular momentum (naturally conserved in spherical symmetry). The first two equations 
above are a rewrite of equation (I3.5P in which dr/dX was expressed from the null condition. 
The third equation is a square root of the null condition. Its integration has to be stopped 
at the ray's turning point where dr/dX = and restarted after that. The sign in the third 
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equation is determined by whether the ray approaches the turning point or goes away from 
it. 

One needs the initial conditions of integration at point F on Fig. 13.11 There the ray 
exits the void forward in time or equivalently enters the void when traced back in time. The 
triangle ABF consists of physical lengths measured by the matter-comoving observer at 
point F for an infinitesimal ray segment. The physical length of the ray segment is dl = \dt\ 
because the locally measured speed of light dl/\dt\ = 1 according to the fundamental axiom 
of General Relativity. The null condition (13. 8ft shows that ABF is a right triangle (simply 
because the coordinates (I2.2p are orthogonal). The angle a is obtained using the end point 
of the previous cheese segment: 

a = &rccos(dir ■ FC) , (3.10) 

— * 

where dir shows the spatial direction of the ray traced back in time and the dot denotes 
scalar product in the global coordinates (13 .4ft . It is easily calculated since the global 
coordinates of F and dir are the end conditions of the previous cheese-segment. One 
can write the equations 

dt = (z + l)d\ 

R'dr 

-. = dt cosa 
Rd9 = dt sina . 

All the functions must be evaluated at F where R = r v a(t) and E(r v ) = 0. After little 
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(3.11) 
(3.12) 
(3.13) 



algebra one gets: 



z{F) = z F 



(3.14) 



L 




r v a(t F ) 
r v a(tp)(zp + 1) sina , 




(zp + 1) sina 



a(t F ) 



where tp is the time at which the ray hits point F. The initial condition used for integration 
is the one for L. 

3.4 Cheese-segments 

The final point of the previous void segment, corresponding to the smallest A when tracing 
back in time, is denoted again with F. For a cheese-segment, the origin of the local 
spatial coordinates can be conveniently chosen to coincide with F since there is no center 
of symmetry analogous to the void center. Rays follow straight lines in the homogeneous 
cheese. The w-axis of the coordinates can be oriented along the ray, u = dir, which sets 
9 = and L = 0. The geodesic equations ( 13.91) simplify to 



dz 

dX 
dt 

dX 

dx 

dX 



z+l 




z + 1 



a(t) ' 



(3.15) 
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where a(t) = (t/to) ' is the EdS scale factor. The ray trajectory in global coordinates is 

x = F + x{\) dir . (3.16) 

— * 

The spatial vector dir gives the direction along the ray traced back in time. The minus 
sign in the third equation above means that x(X) increases when A decreases (going back 
in time), corresponding to going away from the point F in the direction of dir as it should 
be. The initial conditions are the obvious 



z(F) = z F , t(F)=t F , x(F) = 0. (3.17) 

The corresponding analytic solution is: 

u = 1 + ^-^(1 + ^) (3.18) 
otp 

t(A) = t F u 3 / 5 

z (\) = (z F + l)u~ 2 / 5 - 1 

x(\) = i*£(l-„V5). 
dp 

3.5 Tracing luminosity distance 

The area distance by definition is equal to a/ A/Vl s , where Vt s is an infinitesimal solid angle 
subtended by a conical ray at the point source and A is the ray cross-sectional area at 
some distance from the source. If the ray has an opening angle j3 1 at the source, then 
n s = 2tt(1 - cos/3) w 7T/3 2 . For rays with negligible shear, which is the case considered in 
the present paper, the cross section at some distance is a circle of radius /. The angular 
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diameter distance at that position is by definition 




^WS^-' (3 ' i9) 

equal to the area distance. The reciprocity theorem, first derived in 1933 by Etherington 
and popularized by Ellis [29], connects rg (original notation) - the area distance for a 
past directed ray emitted from the Earth observer towards the supernova, to tq (original 
notation) - the area distance for a future directed ray emitted from the supernova towards 
the Earth observer: 

r G = r (l + z), (3.20) 

where z is the redshift between the supernova and the Earth. The luminosity distance to 
the supernova is [29] 

d L =r G (l + z)=r (l + z) 2 . (3.21) 

It is quite surprising that these relations are true in arbitrary spacetimes not only in 
homogeneous ones. 



The angular diameter distance along a past directed ray obeys, 



10|: 



1 d?d a 1 d 2 \/~A 1 , ,,, „ 9 , 

- --R^k v -a 2 , (3.22) 



d A d\ 2 ~ y/A d\ 2 2 

where a is the ray shear, R^ v is the Ricci tensor and k^ 1 = dx^/dX is the null tangent 
vector of the ray. The inverted Einstein equation R^ v = SuGiT^y — g^T/2), where T^ v 
is the energy-momentum tensor and T is its trace, can be used to calculate the righthand 
side of the above equation: 

R^k" = MG{T^k v -^k-k) = 8nG T 00 (k ) 2 . (3.23) 
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The last equality follows from k ■ k = and the fact Tqq = p is the only non-zero component 
for a presureless dust in the comoving coordinates (12. 2p . Finally, k° = dt/dX = z + 1 and 
the angular diameter equation is 

1 d dA -AnGp(X)(z(X) +1) 2 -a 2 . (3.24) 



c^4 dX 2 

The initial conditions at the Earth observer are 



«U(0)=0, ^f°>=-l. (3-25) 



The last equality is a consequence of the relation ddj± ~ where the physical distance 
along the ray close to Earth is dl = \dt\ = \d\\ due to the local speed of light c = \dl/dt\ = 1 
and the normalization choice (13. 2p . The negative sign reflects the choice that A decreases 
tracing the ray back in time, while the distance dj± increases. 

Equation (13.241) has to be integrated along with the geodesic equations on each segment 
of the trajectory. The density p(A) is calculated from the first equation in (12.41) in which 
r = r(A) and t = t(A) are the ones obtained from integrating the geodesic equations. Once 
the angular diameter distance dj± = is known, the luminosity distance is obtained from 
equation (13.211) . 



The shear is generated by inhomogeneous distribution of matter, 

do- , „ / 3M \ L 2 . 

- + 2e e . = A,G P (l-^)- 2 , (3.26) 

where L is the angular momentum integration constant in H3 .91) and 9 e — (1/ ^/~A)d\/~A/dX = 
(l/d^)ddj[/dX is the expansion of the ray bundle. The righthand side is zero in the 
homogeneous cheese and the solution is decaying, a = const d^ . In an inhomogeneous 
void, one can estimate L 2 /R 2 = {d6/dX) 2 R 2 ~ (ii/2r v ) 2 r 2 ~ 1, the bracketed term in 
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(I3.26P is ~ 1 within the void wall and zero outside and the solution is a ~ 47rGpAA with 
the wall thickness AA ~ 5 Mpc. The ratio of the two terms in equation (I3.24p can be 
estimated as cr 2 / \4ttG p) ~ 47rGpAA 2 . The current homogeneous density of the model 
p(to) ~ 10~^ for which AtxG p AA 2 ~ 10^^ and the shear can be safely neglected in ( 13.24]) . 



It plays a significant role only close to mass concentrations like black holes, 



io| 
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4 Distance modulus averaging for a single void 



Discussing a universe with a single void lays the foundation for understanding the 
case of many voids. The Earth observer will always be at the origin of the global spatial 
coordinates. The void center is placed on the X shown on Fig. 14.11 

Y 



The past-directed light rays look unremarkably straight, at least around time £q- This 
is due to the chosen spatial coordinates inside the void (12. 7p . which are approximately 
Newtonian around time £q. The density inhomogeneities are not sufficiently large to cause 
a noticeable light deflection in Newtonian coordinates. 

4.1 Redshift surfaces 

The computed redshift along a given ray is shown on Fig. 14.21 The two "bumps" on 
the plot are where the ray encounters the front and the back of the void. The peculiar 
velocity in the front is towards the Earth observer and the redshift is lower; the velocity 
in the back is pointed away and the redshift is higher. The matter in the EdS cheese has 
a zero peculiar velocity - it does not feel the presence of the void thanks to the matching 
conditions (12. 6p . As a consequence, the redshift outside the void is indistinguishable from 



o 




x 



Figure 4.1: Single void on the X axis 
penetrated by past-directed rays ema- 
nating from the Earth observer at the 
origin O. The minimal cone for averag- 
ing is defined by opening angle (3 max . 
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the one of homogeneous EdS as Fig. 14.21 indicates. The inhomogeneities, through which 
the ray passes, do cause a minute change in the redshift but it is a second order relativistic 



effect 15[ • Primary changes in the redshift result from peculiar velocities of the supernova 
sources inside the void. 

Redshifts in the bumps on Fig. 14.21 correspond to three A positions along the ray. For 
a universe filled with many voids, some rays encounter a double bump when they exit a 
void wall and enter another void nearby. Redshifts in such double-bumps can correspond 
to five A positions. 




Figure 4.2: Solid line - redshift z(X) for 
a ray emitted at 2 ° angle with the X 
axis and penetrating a void centered 
at X = 60 Mpc. Dashed line - redshift 
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ioo ^' MpC z {ty f° r the homogeneous EdS 



All supernovas with the same redshift lie on a spatial redshift surface. A few of those 
are shown on Fig. 14.31 - the actual surfaces are obtained by rotation around the symmetry 
X axis. The numbers that label the surfaces are the corresponding redshifts multiplied 
by 10 4 . Surfaces 70, 75, 203 and 212 consist of three sheets for certain directions - three 
points giving the same redshift along a ray, corresponding to a "bump" on Fig. 14.21 The 
evolution of the triple surface is traced on Fig. 14.31 as the redshift increases to the right, 
two of the sheets of surface 70 merge on the X axis and open up into surface 75 which 
gradually straightens out into surface 128 somewhat before the middle of the void. The 
reverse metamorphosis is seen continuing to the right. Surface 128 folds back into surface 
203 which pinches off into the lense-like and open sheets of surface 212. At some distance 
before and after the void, the redshift surfaces 58 and 235 are just trivial spheres centered 
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X 



Figure 4.3: Redshift surfaces for a void 
centered at 60 Mpc. The bold num- 
bers show z x 10 4 . Some surfaces have 
multiple sheets. The void boundary is 
shown in dashed line. 



at the Earth observer. The surfaces become more symmetric with respect to the middle of 
the void as its distance to the origin increases. 



4.2 Averaging by solid angle and by mass 

The supernova data displays not the luminosity distance d^ itself but a logarithmic measure 
of it, the distance modulus: /i — 5 log^d^/ 10 pc) . The goal of this paper is to compare 
the distance moduli of the inhomogeneous universe and the homogeneous background EdS: 

AvL(z)=iM(z)-iM EdS (z) = 5]ag 10 

where d^(z) is the inhomogeneous luminosity distance (numerically computed by ray trac- 
ing) at redshift z and dj^ dS (z) = 3£q(1 + z — y/1 + z) is the homogeneous one. For small 



df dS (z) 



(4.1) 
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corrections Ad^ = di — (iffS <g; fiff 1 ^ the logarithm can be expanded: 

Ap«(5/lnl0)^, (4.2) 
a L 

showing that to first order the change in the distance modulus is given by the fractional 
change in the luminosity distance. 

In an inhomogeneous universe, d^(z) depends also on the ray direction and so does 
Afi(z) . The correction, averaged over all directions, is obtained by a weighted integration: 

(AM;)> = JM ^ M T (<&) , (4.3) 

where dW(dir) is the probabilistic weight assigned to the infinitesimal element of the spatial 
surface or volume over which the averaging is performed. The element is seen in the 
direction dir from Earth. The probability for observing a supernova emission from that 
element is dP = dW(dtr) / J dW(dir). 

The previous works assume the probability is proportional to the observational solid 
angle on Earth, dW = dQ, and the averaging is over a redshift surface. That is reasonable 
if all light sources are in the cheese where the density is homogeneous and the probability 
for a supernova emission is uniform. For practical calculations, the full 4tt solid angle of the 
Earth observer is broken down into small solid angles Afij. A ray i with direction inside 
AQi is chosen as its representative. The averaging integral is approximated by the sum 

where A/x^ is the distance modulus correction (with respect to the homogeneous EdS) 
evaluated at the points where ray % pierces the sheets of the redshift surface z. One ray will 
usually have several such points and a corresponding number of repetitions of Af^. As a 
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result, the total solid angle will be ^ Af2j > An if the sheets are more than one anywhere, 
which often happens in inhomogeneous universe. The averaging by solid angle becomes 
nonphysical if the redshift surface cuts through a void since a significant probability of 
a supernova emission will be assigned to the void interior which is almost empty. The 
results of this type of averaging will be presented only for reference purposes to show 
that the nonzero corrections are not artifacts of the more physically appropriate averaging 
described below. 

The present paper assumes that the probability for observing a supernova emission from 
a given comoving volume is proportional to the rest mass inside it, dW = dm, in agreement 
with the astronomical views on supernova bias. The rest mass inside the volume specified 
by an observational solid angle on Earth dQ, and sandwiched between redshift surfaces z 
and z + dz is 

E_ pd\ (-2 + 1) 
pdAdl = y^ ,~ , „, dQdz, 4.5 
^ \dz/d\\ 

where the redshift surface area cut off by rays with directions within dVt is dA = d^ dVt 
(from (I3.19p ) and the physical distance between the surfaces is dl = \dt\ = (z + 1) \dX\ = 
(z + 1) \dz\/ \dz/d\\. The sum is over all sheets of the redshift surface, if many. The 
averaging integral is discretized analogously to the previous averaging procedure: 

,a J2i All; Am; £i A ^ Pi d % ( z + ±)AzAQ t / \dz/dX\i 

(Au(z)) = — — — = — (4 6) 

J2i Am i EiPid%(z + l)AzAQi/\dz/dX\i 

The sum runs over all rays and all points where ray % pierces the sheets of redshift surface 
z. 



30 



4.3 Maximal average correction in the minimal cone: numerical 
results 

The real universe does not contain vast homogeneous regions in between the voids. Such 
regions have A/i = and will damp the correction coming from the inhomogeneous void 
when averaging is performed. The largest possible average (A/a(z)) is obtained when the 
amount of cheese is minimal. With that in mind, the averaging for a single void is restricted 
to rays inside the minimal cone with an opening angle Pmax on Fig. 14.11 The minimal 
cone contains the void and a minimal amount of cheese. The obtained average corresponds 
to the best case scenario when all the voids in that range are exactly at the same radial 
distance from the Earth observer, forming a spherical shell, and their average corrections 
add up constructively to the total average. 

The void is centered at distance d = 60 Mpc on the X axis. Due to the axial symmetry, 
it is sufficient to shoot rays only in the upper half of the XY plane and the full picture is 
obtained by rotation around the X axis. The angle from zero to f3 ma x — arcsin(r w / d) — 30 
was divided amongst 3000 rays separated by A/3 = 0.01 °. Ray i represents the solid angle 
Afij = 27i(cos(Pi — A/5) — cos((3i)) between the cones with opening angles /3j — A/3 and /%. 
The redshifts for averaging were 500 values equally spaced on the interval 0.001 < z < 0.025. 
For each redshift, the points where the rays intersect the redshift surface were found and 
the average correction (Afi(z)) was calculated according to (14.41) and (14.61) . The results are 
shown on Fig. 14.41 

Although the solid angle average is not physically sensible for sources inside the void, it 
has the same general behavior as the mass average. This shows the results of averaging are 
qualitatively independent of the chosen averaging procedure. The redshifts with positive 
average corrections correspond to surfaces before 128 on Fig. 14.31 which are shifted forward 
with respect to the EdS surfaces, thus having a bigger d^. Surface 70 (redshift 0.0070) 
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^ ^_ Figure 4.4: Correction (Afi(z)} aver- 

°- 025 aged by solid angle (dashed) and by 

mass (solid) for a void centered at 60 

Mpc. 



is the one with the maximal mass-averaged correction on Fig. 14.41 Surfaces after 128 
have negative corrections since they are shifted backwards with respect to the EdS redshift 
spheres. The smaller magnitude of the negative corrections compared to the positive ones is 
caused by the bigger amount of cheese in the minimal cone at higher redshifts which damps 
the average towards zero. The positive and negative corrections become more symmetric 
with each other for voids at larger distances from the origin - they are completely symmetric 
in the limit of an infinite distance. 

The small spikes appearing on Fig. 14.41 and later figures are errors from the numerical 
discretization f !4.6p of the integrals in (14. 3p . Notice they appear at redshifts in the proximity 
of surfaces z = 0.0084 and z = 0.0186 on Fig. I4.3l which have a segment pointing straight to 
the observer. It is hard to integrate over such a segment with finite number of rays in (14. 6p . 
Increasing the number of rays decreases the oscillations but becomes time consuming. The 
reader should simply ignore the oscillations. 

The observer may want to bin the data and find the average correction, [A/x], over a 
redshift interval. Naively, one can sum up all the 500 corrections on Fig. 14.41 and divide by 



their number to obtain [A/i] 



naive 



0.014 for the angle-averaged and [A/i] 



naive 



0.010 for 



the mass averaged correction over the whole plotted interval. More appropriately, one has 
to weigh each (Afi(z)) by the total weight of its redshift surface z, shown on Fig. 14.51 and 
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Figure 4.5: Total solid angle for each 
redshift surface used for averaging. 
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Figure 4.6: Total rest mass for each 
redshift surface used for averaging. 



Fig. 14.61 The weights on Fig. 14.61 are increasing because there is more mass in the minimal 
cone at larger redshift. The two peaks in the front and back of the void are again due to 
overcounting caused by multi-sheet redshift surfaces there. The weighted averages over the 
redshift interval shown on Fig. 14.41 are [A/i] = 0.016 for angle averaging and [A/i] = —0.006 
for mass averaging. The interval averages [A/i] are significantly smaller in magnitude than 
the single-redshift averages (Afi(z)) on Fig. 14.41 That results from the tendency of the 
positive and negative corrections {Afi(z)) from the front and back of the void to cancel out. 
The cancellation is imperfect close to the Earth observer since the redshift surfaces are not 
exactly symmetric with respect to the middle of the void. For voids at larger distances, the 
symmetry improves and so is the cancellation. A perfect symmetry is achieved when the 
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void is at infinite distance from the origin and the rays run parallel to the X axis. Then 
the interval average would be exactly zero. If one wants to avoid the corrections shown on 
Fig. 14.41 the data must be averaged over intervals bigger than a single void. 

The interval averages calculated above implicitly assume that the supernova detection 
efficiency is constant throughout the redshift interval. That is usually not true in obser- 
vational astronomy since supernovas at higher redshift are harder to detect and less in 
number. The detection efficiency would be especially non constant for big redshift inter- 
vals/big voids. That would bias the interval average [A/j] towards positive values. 

4.4 Maximal average kinematic correction vs. photon number 
conservation 

The argument [12j based on the gravitational lensing conservation of the total photon 
flux states that the luminosity distance, averaged over all directions in an inhomogeneous 
universe, is the same as in the homogeneous counterpart. The inhomogeneities induce 
focusing/defocusing of light but the argument ascertains these are random along different 
lines of sight and cancel out in the average. In reality, they are correlated close to the 
observer and the average does not have to be zero. 



The argument in 12| is implicitly staged in the Newtonian coordinates (12. 8ft . It assumes 
that the redshift surface in the inhomogeneous universe remains spherical and at the same 
coordinate distance from Earth as in the homogeneous counterpart. The analysis in the 
previous sections shows this is not the case because the matter peculiar motion can shift 
and bend the redshift surface as shown on Fig. 14.31 If the redshift surface shifts along a 
ray by a comoving distance As in Newtonian coordinates, the kinematic correction to the 
luminosity distance is 

Ad L « (1 + z)As , (4.7) 
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There are additional corrections to d,£ due to weak lensing or non-linear effects like Rees- 
Sciama. Neglecting them leads to a very good fit to the numerically calculated maximal 
average correction at low redshifts which proves the kinematic correction ( 14.7ft dominates 
there. The other effects become apparent at high redshifts where the kinematic correction 
has already decayed and for big voids into deep nonlinear regime or with a decaying mode. 
Equation ( 14. 2 p shows the distance modulus correction Ajx is proportional to the fractional 
change in the luminosity distance Ad^/dj^ ~ As/s (s is the comoving distance to the 
redshift surface along the ray). It will decrease as A/i ~ 1/s, because the kinematic shift 
As is bounded from above by the void comoving radius. For distances much bigger than 
a void radius, the photon conservation argument will apply but only in the approximate 
sense that with the distance s increasing, the redshift surface is getting closer to spherical 
and the fractional corrections are getting smaller. 

From the linear perturbation formula ( 12.141) . the shift of the maximal correction surface 
70 on Fig. I4.3l is As max oc v max /H ss ar-\8\/3 = r_|5|/3(l + z), where r_ = 25Mpc is the 
radius of the underdense void interior. For models with a linear growing perturbation in 
an EdS background \8\ = \5q\cl, where \Sq\ = 0.9 is the underdensity in the void interior at 
present time. The change in the luminosity distance is Ad^ rj {l + z)As max oc r_|<5o|/3(l + 
z). The corresponding maximal average correction for a given redshift is obtained from 

(Afi(z)) max = 5 log 10 

where (z) = 3£o(l + z ~~ + z ) is the homogeneous luminosity distance. The pa- 
rameter rj represents the ratio of the maximal average correction to the naive maximal 
correction corresponding to the maximal peculiar velocity. The value of rj is not derivable 
from linear perturbation theory but requires a numerical calculation. It depends on the 
details of the model (peculiar velocities) and the specifics of the chosen averaging proce- 



?7|5o|r- 



3(1 + z)df dS (; 



(4.8) 
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dure (mass distribution). It was found numerically that r\ does not depend on redshift for 
the models with a pure growing mode considered here. Models that contain a mixture of 
growing and decaying mode have a more complicated dynamics; the redshift surface shapes 
have a strong time dependence and n ^ const. 

For small corrections, the logarithm in (14. 8 j) can be expanded, showing the correction 
decays like r_/s or in other words oc 1/N V where N v = s/(2r v ) is the number of void 
diameters that equal distance s. 
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Figure 4.7: 
Solid curves 
mass-averaged 

(fi- fi EdS ^ for 
a void centered at 
X =(60, 120, 180, 
240, 300, 360, 420, 
480, 540, 600) Mpc. 
Dashed curves 
maximal average 
correction (14.81) for 
r] = 0.20 (upper) and 
r\ = —0.18 (lower). 

Dotted curve shows 
^LCDM _ ^EdS 



Fig. 14.71 is a composite plot showing the distance modulus correction averaged by 
mass in the minimal cone for a void placed at various distances from the Earth observer. 
The void positions do not intersect each other but the corrections overlap since a void 
can influence a nearby redshift surface that is outside - surfaces 70 and 212 on Fig. 14.31 
gained extra sheets because of the void. The positive and negative corrections get more 
symmetric at bigger distances from the observer. The dashed curves on Fig. 14.71 show the 
maximal average correction estimated by (14.81) . The upper curve, enveloping the maximal 
positive corrections, corresponds to r\ — +0.20 and the lower curve, describing the negative 
corrections, is for rj = —0.18. At small redshifts, the voids enter nonlinear regime and the 
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correction is higher than the linear estimate (14. 8p . The dotted curve on Fig. 14. 71 corresponds 
to the standard LCDM with Q m = 0.3 and Q\ = 0.7 - obviously the average correction 
cannot substitute for dark energy even in the limited redshift interval considered. 

The simple estimate (14. 8 p of the ideal scenario maximal correction does not take into 
account the presence of an excessive amount of cheese inbetween the voids nor the fact that 
the corrections from randomized voids often add destructively as will be seen later. The 
first effect is an artifact of the Swiss-cheese model - the real universe does not contain vast 
homogeneous areas between the voids. Both effects will damp the averaged corrections of 
the model and they will go to zero with redshift faster than shown on Fig. 14.71 
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5 Cumulative distance modulus correction by aligned 
voids 
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Figure 5.1: Correction Afi for a ray 
at 0.02 degrees with the X axis pen- 
etrating a string of 60 voids centered 
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at X = 2r v i with i 
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The cumulative distance modulus correction for a ray at 0.02 degrees angle piercing a 
string of 60 voids (r v = 30) aligned along the X axis is shown on Fig. 15.11 The voids 
touch but do not intersect. Each dot shows the correction right after the ray exits a void. 
The correction inside the voids is not shown - it is multivalued for some redshifts due to 
multi-sheet redshift surfaces. Including the ray shear in the luminocity tracing equation 
does not influence the result. A correction of a similar magnitude at large redshifts was 

n 

obtained in [9] but for much larger and significantly nonlinear (density contrast in the void 
wall ~ 28) voids of radius 350 Mpc. The effect was ascribed to the cumulative gravitational 

n 

lensing from the underdense void interiors [ll| and shown to go to zero when averaged over 
randomized impact parameters. Fig. 15.11 demonstrates that a large correction is possible 
for small voids too. The cumulative correction at a given redshift does not depend strongly 
on the size or the nonlinearity of the voids but only on how empty they are along the ray. 
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One of the main questions addressed in this paper is whether the significant cumulative 
correction at big redshifts persists in the distance modulus when it is averaged over all 
directions in a lattice of voids. The answer is no, as shown in the next section, and is 



The 



obtained without any approximations or assumptions unlike the previous studies 
more surprising fact is that randomization of the voids is not necessary - even regular 
void lattices like the simple cubic have a vanishing average correction at high redshifts. 
Of course regular lattices will still show large cumulative corrections in certain directions, 



as 30| demonstrates, but the required void alignment is extremely improbable in the real 



universe. 
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6 Averaging in lattices of r v = 30 Mpc voids 



The real universe is permeated by a lattice of voids. Its effect on the averaged distance 
modulus is calculated by shooting past-directed rays from the Earth observer in directions 
covering the full solid angle. The luminosity distances at given redshifts are averaged over 
all directions via (14. 6ft . 

A lattice of 1691 voids was obtained by randomly generating void center coordinates 
(x, y, z) on the interval —390 < x,y, z < 390 Mpc. The voids can touch but not intersect 
each other. They are not allowed to contain the Earth observer or get closer than 1.5 r v to 
her since that would lead to too big (A/x) corrections. The closest void center is at distance 
d = 46.9 Mpc from Earth. The amount of unwanted cheese between the voids depends 
on the packing efficiency of the lattice, defined as the total volume fraction that is inside 
voids. The maximal packing efficiency of identical spheres is n/y/18 ~ 0.74 achieved in the 
face-centered cubic and the hexagonal close-packed lattices. The generated random lattice 
has a packing efficiency of 0.35 implying too much cheese. Unfortunately, higher packing 
cannot be achieved by the simple random process described here. 

For comparison, the distance modulus correction was also averaged in a simple cubic 
void lattice with a cube cell size of 2r v and a packing efficiency of n/Q w 0.52. In this case, 
one can average only over the first octant, the others being identical, thus reducing the 
computational work 8 times. The closest void in that lattice is at distance d = \/3r v m 51.96 
Mpc from the Earth observer. 

For the random void lattice, the total observational solid angle of 4 7r was divided 
approximately equally between 208024 directions/rays. The resolution in the observational 
polar angles was d6 = d<fi = 2/257 i.e. the rays could resolve a detail of a comoving 
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size 2 Mpc at a comoving distance 257 Mpc (corresponding to redshift z = 0.063). The 
luminosity distance was traced along each ray and its values corresponding to 200 redshift 
points equally spaced on the interval < z < 0.08 were extracted. Those values were 
used in (14. 6[) to calculate the average for each redshift point. The resulting mass- averaged 
correction is shown on Fig. 16.11 For the simple cubic lattice, the 47r/8 solid angle of 
the first octant was divided between 83599 directions/rays. The chosen resolution was 
d9 = dcj) = 1.5/345 i.e. the rays could resolve a detail of size 1.5 Mpc at comoving distance 
345 Mpc (corresponding to redshift z = 0.086). The redshift points were the same as for 
the random lattice and Fig. 16.21 shows the obtained average correction. The thin solid line 
shows the standard deviation of the corrections at each redshift a[Afi](z) = ^Jvar [A/x] with 
the variance var [A/i] = (A/i 2 ) - (A/i) 2 . The averages (...) over directions were calculated 
with the mass weights dW = dm. 



Figure 6.1: Thick solid curve - mass- 
averaged (A/x) for the random void lat- 
tice. Dashed curve - maximal average 
correction for rj = 0.20. Thin solid 
000 002 004 006 008 curve - standard deviation of A/i at 

z each redshift. 

The average corrections on Fig. I6.1l and Fig. [672] are damped below the maximal average 
correction (14.81) given by the dashed curve. The first reason for that is that the two lattices 
have packing efficiencies lower than the one in the minimal cone, which approaches ~ 2/3 at 
distances bigger than several radii. The bigger amount of cheese with zero peculiar velocity 
in the lattices compared to the minimal cone damps the average correction towards zero. 

The much more serious reason for the damping is that whenever a void front and 

another void back are at the same distance from the Earth, their corrections will tend to 
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08 Figure 6.2: Same as the previous figure 
but for the simple cubic lattice. 



Figure 6.3: Distribution of voids in radial 
distance (Mpc) for the random lattice. 



Figure 6.4: Same as the previous figure but 
for the simple cubic lattice. 



cancel out in the total average. The astonishing efficiency of this cancellation damping, 
especially in the random lattice, at redshifts even as small as a void diameter is puzzling. 
This can be understood visually by looking at Fig. 16.31 and Fig. 16.41 which show the void 
distribution in radial distance for the random and simple cubic lattices. Obviously, the 
front-back cancellation is more severe in the random lattice and at higher redshifts since 
there is more volume available there and correspondingly more voids per redshift interval. 
Quite surprisingly, Fig. 16.21 demonstrates that the cancellation occurs very vigorously also 
in a regular void lattice and even there the cumulative corrections from the previous section 
do not survive in the directional average. 

The voids in the real universe are not exactly spherical which allows for a tighter packing 
and more correlated radial void positions close to the observer than in the latices considered 
here. The expected damping would be between that of the random and the cubic lattice 
which in this way serve as lower and upper estimators of the average correction. Further 
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discussion of the numerical results is given in the summary and conclusions section. 
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7 Averaging in a lattice of r v = 300 Mpc voids 



A natural question to ask is whether bigger voids would produce a proportionally bigger 
average correction than the one obtained so far. 

The current astronomical data suggests that there are structures of size significantly 



larger than 30 Mpc. Galaxy count studies 



3 if ] imply the Southern local void may have a 



radius of 200 Mpc and a density contrast of —0.25. The size and underdensity of the largest 
voids is restricted by the imprint they would leave on the CMB passing through them. A 



recent paper [32j found such an imprint by averaging the CMB temperature around voids 
at a mean redshift z v = 0.53. Initially, the signal was attributed to the late-time integrated 
Sachs- Wolfe effect. Since a significant linear Sachs- Wolfe effect is absent in a flat matter- 



only universe (see references in 



33|). the signal was presented as a direct manifestation of 



34j showed that the measured signal is 



dark energy. However, two subsequent papers [33|, 
orders of magnitude too large to be explained with the Sachs- Wolfe effect in LCDM. 

The matter density of the voids in this section is chosen to reproduce roughly the average 



imprint on CMB, shown as a dashed line on Fig. 9 in 33]. The necessary matter density 



shown on Fig. 17.11 is 



A\ + A 2 tanh[0.2(r - 200)] - A 3 tanh[0.3(r - 290)] , r < 300 
p(r,t ) = p(t Q ) { ' (7.1) 

1 , r > 300 



with the coefficients (A 1: A 2 , A 3 ) = (0.764715, 0.349972, 0.115257), determined by the 
mass-compensation condition (12. 6p and density continuity. A more precise match to the 



curve in 33| would require a density function with more parameters and is not the goal 
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of the present paper. The bang time chosen for the model is £g(r) = 0, implying a pure 
growing mode. The homogeneous "cheese" between the voids is the same Einstein-de Sitter 
metric that was used for the small voids. 

p (r, t Q ) I p (t ) 
1.2- 



Figure 7.1: Density ratio p/p at time 
50 ioo 150 200 250 300 350 *' MP ° *0 for a void of radius r v = 300 Mpc. 



a, 



3 




4 6 8 10 

angle (degrees) 



Figure 7.2: Change in the CMB tem- 
perature (pK) produced by a void of 
14 radius r v = 300 Mpc centered at red- 
shift z v = 0.53. 



Fig. 17.21 shows the effect of such a void, centered at redshift z v = 0.53 (1633 Mpc in 
EdS), on the CMB temperature. It is calculated by shooting rays as shown on Fig. 14.11 

and propagating them to some redshift zp and time tp outside the void. The redshift is 

3/2 

propagated further to the time of last scattering tpg = tQ a^g using the analytic solution 
( I3.18P : zls = (zp + IX^fAlsO 2 ^ 3 — 1- The fractional change in the CMB temperature 
is AT/T = —Az/z = —{zpg — z ps)/ z ls> wnere Z LS = 1099 i s the unperturbed redshift 
corresponding to apg = 1/1100. 



The effect shown on Fig. 17.21 and in 



33] is orders of magnitudes smaller than the one 
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calculated in [30|. The large AT/T in [30|] is due to: (1) a cumulative effect along special 
directions from aligned voids in a regular lattice, (2) a significant gravitational blueshift 
of the rays starting from inside voids at the initial time of the model when the density 
contrast in the void interior is practically S{t^ t ) = — 1, (3) possibly effects produced 
by the significant non-linearity of the chosen void model jj] having a density contrast in 
the void wall 5 (to) — 28 In the real universe, special alignment of more than 2 - 

3 voids to produce a cumulative effect on CMB has a vanishing probability. The CMB 
radiation starts its journey towards us when the dark matter density contrast is of the 
order \5\ ~ 1CP^ producing a gravitational blueshift much smaller than the one found in 



I ■ 

301 ] . These considerations contradict the claim in [30| that voids of radius larger than 



y 



the large-scale part of the CMB 



22| in the SDSS data. 



35 Mpc are excluded because they modify significant 
temperature spectrum. Such voids are actually found 

The cumulative correction \i — ii E ^S a i ori g a ra y a t 0.02 degrees with the X axis, 
penetrating a string of 6 non- intersecting voids centered at X = r v (2i — 1), % = 1...6 is 
A/i = 0.11 at redshift z = 1.97. This is about twice as low as the correction for the small 
voids (A/x = 0.27), see Fig. 15.11 The density contrast in the interior of the bigger voids 
is also about twice as low as for the small r v = 30 Mpc voids which indicates that indeed 
the cumulative effect at a given redshift does not depends on the void size but on how 
empty the voids are along the ray. Analogously to the case of small voids, the cumulative 
correction is completely destroyed in the averaging over all directions, even for the regular 
simple cubic lattice. 

Fig. 17.31 shows the correction averaged in the minimal cone for several positions of a 
single void. It is counterintuitive that the correction at distance r v is not bigger than the 
corresponding one for the small void. Its magnitude is determined by the ratio Adj^/vy oc 
v max/ r -- The explanation is that although the bigger voids have proportionally larger 
peculiar velocities, the linear formula ( 12. 14ft shows that the ratio v ma x/r- is independent 
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Figure 7.3: Solid 

averaged ( jj, 



curves - mass- 
for a void of 



radius r v = 300 Mpc centered at 
X =(400, 800, 1200, 1600, 2000, 2400) 
Mpc. Dashed curves - maximal av- 
erage correction for 77 = 0.22 (upper) 
and 77 = —0.20 (lower). Dotted curve 
shows n LCDM - [i EdS . 



Figure 7.4: Thick solid curve - mass- 
averaged (A/i) for the simple cubic 
void lattice. Dashed curve - maximal 
average correction for 77 = 0.22 and 
r v = 300 Mpc. Thin solid curve - stan- 
dard deviation of A/x at each redshift. 



of the characteristic void size r_. The dashed curves on Fig. 17.31 are the maximal average 
correction estimate (14. 8 j) for r_ = 200 Mpc and \8q\ = 0.47. The linear perturbation theory 
estimate (14. 8p envelops the maximal correction solid curves very well except at low redshifts 
where the voids become mildly nonlinear. It is a little surprising that approximately the 
same values of rj from Fig. 14.71 work on Fig. 17.31 for an entirely different void model. 
That suggest some kind of universality of those values for the mass averaging procedure, 
at least for mass-compensated voids that have a pure growing perturbation mode, are at 
most mildly nonlinear, and have approximately constant densities in the void interior and 
walls. 

Fig. 17.41 shows the correction averaged over all directions in the simple cubic lattice. 
The random lattice was not used since it underestimates the average by having too much 
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cancellation damping as was discussed in the case of small voids. To calculate the average 
(A/i), the simple cubic lattice used for the r v = 30 Mpc model was rescaled to the new void 
radius r v = 300 Mpc by multiplying the void center coordinates by 10. The 4 tt/8 solid angle 
of the first octant was divided between 51126 rays with a resolution d6 = d<fi = 10/1800 
at distance 1800 Mpc (redshift 0.60). The redshift points for averaging were 200 and 
equidistant on the interval < z < 0.8. The average correction (thick solid line) is damped 
with respect to the maximal average correction (14. 8 j) (dashed line) due to cancellation 
between fronts and backs of different voids, analogously to the case of small voids. Further 
discussion of the numerical results continues in the summary and conclusions section. 
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8 Summary and conclusions 



The paper studied how the distance modulus, averaged over all lines of sight in an 
inhomogeneous universe, differs from its homogeneous counterpart. The inhomogeneities 
were represented by random and regular lattices of mass-compensated voids in Swiss-cheese 
models. Two void layouts were considered with a small (30 Mpc) and a big (300 Mpc) 



22j and the second deduced from its imprint on 



radius, the first observed in the SDSS data 
CMB 32| . The Earth observer was put in the cheese but the conclusions will not change 
significantly for an observer inside one of the voids. For the first time, the averaging 
of the distance modulus was widened to include the supernovas inside voids. That was 
made possible by assuming that the probability of supernova emission from a comoving 



11] . the average 



volume is proportional to the rest mass in it. Unlike previous studies 
correction to the distance modulus was calculated by exact numerical ray tracing without 
using assumptions about the ray impact paramaters or perturbation theory approximations, 
and the result includes all linear and non-linear effects. 

The distance modulus correction, due to gravitational lensin g, t hat accumulates along 



10 



11] was calculated. The 



a ray crossing diametrically a sequence of aligned voids 
correction increases linearly and becomes quite significant at high redshifts, see Fig. 15.11 
It is roughly proportional to the density contrast in the void interior but does not depend 
on the void radius. It was demonstrated that the distance modulus, averaged over all 
directions in a lattice of voids, does not show a cumulative correction at high redshifts. 
That is true not only in a random lattice (see Fig. 16. ip but also in a regular simple cubic 
one (see Fig. 16.21 and Fig. \7A\i implying that void randomization is not necessary for the 
destruction of the cumulative correction when averaging, contrary to popular belief. A 

49 



arg e cumulative correction is still observed along special directions in a regular void lattice 
30J, but the probability for the required void alignment is vanishing in the real universe. 

At low redshifts z, perturbation theory predicts that the radial peculiar velocities v r are 
capable of producing significant fluctuations in the luminosity distance Ad^(z) / 'd^ ~ —v r /z 
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36[ (v r measured in speed of light units). That leads to a scatter in the local Hub- 
ble diagram AH/H = —Ad^/d^ and, by ( 14. 2p . to a non- vanishing average correction 
(Afi(z)) ~ —2.17 (AH/H) ss — 2.17 (v r (z)) j z, where the so called peculiar velocity 
monopole (v r (z)) is the radial peculiar velocity at redshift z, averaged over all directions. 
A measure of the fluctuation magnitude is its average within a sphere of radius R around 
the observer, [AH/H]r, which varies with location. Its cosmic variance over all possi- 
ble observers/locations has been calculated previously in several cosmological models by 



N-body simulations 



trum 
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37l | and by linear perturbation theory using the matter power spec- 



391 ] . Such variances are employed to predict confidence intervals for the value of 



[AH/ H]r measured by a random observer. Since the obtained cosmic variance of [AH/ H]r 



decreases with R 
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39| | , with the cosmic mean being zero, the likelihood to measure large 
average fluctuations decreases with redshift, as expected. 

The cosmic variance of the sphere- average [AH/H]n over all observers cannot be used 
to estimate the average over all directions (AH(z)/H) in the redshift bin z, experienced 
by a particular observer. The results in the present paper allow to calculate (AH(z)/H) 
by studying such an observer in a simple model of our local neighbourhood. The averages 
and variances so obtained are not cosmic, thus are more applicable to our astronomical 
observations. They cannot be calculated from a given matter power spectrum but require 
a concrete numerical realization of an inhomogeneous universe. To the best knowledge 
of the author, this is the first time when the non-vanishing direction-averaged correction 
(Afi(z)) is calculated as a function of redshift for universes containing voids. The upper 
bound (14. 81) for (Afi(z)) was motivated using linear perturbation theory and was found to 
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agree with the numerical calculations on Fig. 14.71 and Fig. 17.31 The parameter r/ depends 
on the particular form of inhomogeneities and the choice of averaging procedure. The 
numerically obtained value rj ~ 0.20 indicates that the maximal average correction is 20% 
of the naive correction corresponding to the maximal peculiar velocity. This rj turned out 
to be a universal constant, approximately independent of the void size, for voids which 
have approximately constant densities in their interior and walls that are not in a deep 
nonlinear regime. Due to void randomization, it is expected that the average correction 
will be suppressed below the maximal 20% and will tend to zero. The efficiency of that 
process and at what redshift it sets in has not been studied before. The calculations 
revealed that the average correction within void lattices is indeed severely damped below 
the predicted maximal average correction due to cancelations between the fronts and backs 
of different voids. As figures 16. 1\ \6.2\ and 17.41 demonstrate, the cancelation is surprisingly 
efficient at low redshifts even in regular lattices and the average correction drops below 
0.01 mag after a single void diameter. Nevertheless, the average correction is not zero close 
to the observer, indicating that the implicit assumptions of the photon flux conservation 



argument stated in 12] do not apply at low redshifts. 

The nonzero (A//(z)) is observable, provided there are enough supernovas N z in redshift 
bin z to reduce the sampling statistical error of the average: a[ (A(j,(z)) } = a[A/i(z)]/ 'y/Nz, 
where a[Afi(z)} is shown as a thin solid curve on figures ISTTl 16.21 and 17.41 An underdense 



"Hubble bubble" around us with a Hubble parameter s 
(correspondingly (Au(z)) < 0) has been argued for in 



ightly higher than the global one 
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not confirmed in 
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411 ] . However, that claim was 



43j and was shown to depend on the way the supernova magnitudes 



were extracted from the photometric data 



44|. Local bubble or not, fluctuations in the 



distance modulus binned average begin to appear in the newest data sets, Fig. 20 in [451 ]. 



and Fig. 7 in 



46] 



Significant efforts in contemporary cosmology are aimed at measuring a possible time 
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evolution in the dark energy equation of state. That would require fixing the Hubble 
constant at low redshifts to a 1% precision [dj], 47]. Ongoing and future low redshift 
surveys will boost the number of the observed nearby supernovas sufficiently to bring the 
sampling error of the average a[(AH / H)] « er[ (Afi(z)) ] / 2.17 below 1%. In contrast, the 
directional average itself (AH/H) w — (A/i(z)) / 2.17, generated by the coherent peculiar 
motion, is not influenced by the supernova number and will degrade significantly the errors 
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48| . Consequently, the peculiar velocities 



in measuring the dark energy time dependence 
need to be subtracted from the low redshift Hubble diagram. 

Several methods are utilized for that. The Local Group (LG) of galaxies moves as a 



whole with respect to CMB at a speed v^q, the so called peculiar velocity dipole 



50, 



5 if ] . The most common way to correct the Hubble diagram for that motion is to adjust 



the observed redshifts to an observer riding the LG barycenter. The luminosity distance 



fluctuations seen in that frame are Ad^/d^ ~ — (v r — v^q)/z [35|] and the velocity dipole 
will cancel out since the coherent velocity component of a nearby galaxy with respect to 
CMB is v^ 1 ~ vlq. Not correcting the CMB redshifts to the LG frame is permissible for 
a large supernova sample that covers densely and uniformly all directions: the coherent 
dipole fluctuation seen in the CMB frame is Ad^/d^ —v r /z ~ ~^LG cos {@) l z an d 
that expression vanishes when averaged over the full solid angle. A further refinement 
is to correct the LG redshifts for infall velocities in the LG frame caused by the nearby 
superclusters represented by a crude linear multi-atractor model. That was done in the 



Hubble Key Project 49( which measured the Hubble constant to a 9 % precision, but it 



was not very efficient at reducing the scatter in the low redshift Hubble diagram as seen 



on Fig. 4 of [52|. Accordingly, the philosophy of the Hubble Key Project was to extract 
the Hubble constant mainly from secondary distance indicators at high redshifts. A finer 
method to correct for peculiar velocities is to calculate them from the observed galaxy 
distribution utilizing linear perturbation theory and choosing a galaxy-dark matter biasing 
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parameter that minimizes the scatter on the Hubble diagram. That technique does reduce 
the scatter [53] but using it to infer the Hubble parameter from low redshift data produces 
a controversely high value of Hq = 85 km s -1 Mpc~^ j&4] . To avoid unreliable velocity 



45|, 
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55| simply chose to include only objects above 



corrections, many recent surveys 
a certain redshift, z > z m i n , where the peculiar velocities are considered insignificant 
compared to the Hubble flow. The Hubble constant and the dark energy equation of 
state parameter w inferred from the low redshift data are sensitive to the choice of z m j„ 



411, 145]. As discussed in [45], there is a wide variation in the values selected for z m { n in 



the literature, ranging from 0.01 55] up to 0.023 47J, which indicates that there is not a 
universally accepted prescription. 

Figures IBTTl 16.21 and 17.41 of the present paper allow to select z m { n readily. The required 
AH/H = 0.01 precision corresponds to Ang 0a i = 0.0217. The sampling error of the average 
is a[An(z)]/ y/Nz, where cr[A[j,(z)] is given by the thin solid curves on the plots and N z is 
the number of objects/supernovas in redshift bin z. Provided there is sufficient statistics, 
the sampling error will fall below Afig 0a i. In that case the limiting value z m i n is determined 
by the average (Afi(z)) itself - a natural choice is the redshift at which the maximal aver- 
age correction (14. 8 j) (the dashed curves) drops below Afi goa i or the redshift corresponding 
to a void diameter (where (Afi(z)) ~ 0.01 due to front-back void cancellation), whichever 
is lower, in co_, estates based on the sphere-average { AH,H ]h are n,ore pes- 
imistic: its cosmic variance in the LCDM model drops below 1% at z m i n = 0.05 [39] (linear 
estimate). If our neighbourhood contains voids of the size observed on the CMB imprint 



321 ]. the required z m i n would be much larger, as Fig. 17.41 indicates. 

Excluding the objects below z m ^ n reduces the size of the sample and proportionally 



increases the statistical error j48] . The present paper suggests it is possible to preserve the 
low redshift data, instead of throwing it away, by collapsing it into an interval average [A/j] 
over a redshift interval containing a void diameter. These averages turn out very close to 
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zero: [A/i] = 0.003 (z = . . . 0.025 on Fig. EZTJ), [A/i] = 0.007 (z = . . . 0.025 on Fig. 
16. 2p . [A/x] = 0.002 (z = . . . 0.25 on 17. 4p . That stems from the the cancellation between 
positive and negative lobes in (Afi(z)) and the mass averaging giving higher weights to 
higher redshifts where (A/i(z)) is closer to zero. A vanishing [A/i] ~ means [/i] ~ . 

Wi (fi(zi)) ~Yl W i $ dS (H , Zi), (8.1) 

i i 

where Wi is the mass weight assigned to redshift bin Zi, (/i(<2i)) is the distance modulus 
in that bin averaged over all directions, and the sum is over all the bins in the redshift 
interval. The lefthand side of the formula is calculated from the observational data. The 
weight Wi is proportional to the rest mass inside the redshift bin which can be estimated 
as the corresponding mass in the homogeneous background model. Unlike a traditional 
Hubble diagram that weighs all redshift bins equally, here Wi oc zj Azi at low redshifts. 
The righthand side of (18. II) depends solely on the Hubble parameter Hq at low redshifts 
and on additional cosmological parameters of the background model at higher redshifts. It 
allows one to extract the value of Hq from the low redshift data. This method will become 
possible for future surveys that: (1) have a dense full-sky coverage to calculate the average 
correction over all directions (A/i(z)); and (2) contain a large number of objects in each 
redshift bin so that (A/i(z)) is readily observable, not swamped by the sampling error. 

If the observer is inside a void, the average correction (Afi(z)) starts with a negative 
lobe unlike the so-far considered case of an outside observer. A simple illustration of that 
is an observer at the center of a void, which will measure A/i(z) = /i(z) - fi EdS {z) < for 
supernovas inside the void since the matter there is expanding faster than the background 
and it takes a smaller distance to achieve the same redshift. The interval average [A/i] 
over a void diameter would still be very small since the mass averaging scheme adopted 
in this paper ascribes lower weights to low redshifts. Additional studies are needed to 
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investigate the behavior of [A//] in the presence of superclusters - the other type of major 
inhomogeneities encountered in the universe. A plausible guess is that [A/x] vanishes on a 
redshift interval encompassing a supercluster, justifying the validity of ( 18. ip in that case as 
well. 
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